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Abstract. We have observed the dark cloud Barnard 68 with SCUBA at 850 fj.m and with 
SIMBA at 1.2 mm. The submillimetre and millim etre dust emission correlate well with 
the extinction map of Alves. Lada. & Lada < 200 iL The Av/850 yum correlation is clearly 



not linear and suggests lower temperatures for the dust in the inner core of the cloud. 
Assuming a model for the temperature gradient, we derive the cloud-averaged dust emis- 
sivities (normalised to the V-Band extinction efficiency) at 850 yum and 1.2 mm. We find 
i^sso/jm/xw = 4.0 ± 1.0- 10"'' and Ki 2mmlK\ = 9.0 + 3.0- 10"'. These values are compared with 
other determinations in this wavelength regime and with expectations for models of diffuse 
dust and grain growth in dense clouds. 
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1. Introduction 

Despite being a fundamental parameter in Far Infrared, Submillimetre and Millimetre astron- 
omy, few measurements of the dust emissivity' are available (for a review, see 

Send offprint requests to: S. Bianchi, sbianchi@arcetri.astro.it 
' The emissivity proper is the emission [absorption] cross section normalized to the geometrical cross 

section of a dust grain. Other authors prefer the cross section normalized to the mass of the grain. Both 

quantities are the same when normaUsed to the analogous quantity for V-band extinction. 



Alton et all200C ; 
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Fig.l. SCUBA map at 850 jum (left) and SIMBA map at 1.2 mm (right) of Barnard 68, with 
superimposed Ay contours. The Ay and SCUBA maps have been smoothed to match the SIMBA 
resolution (FWHM=24"; the beamsize is shown on the right image). All images have been re- 
sampled to a pixel size of 12". The field of view is 5'x5'. Ay contours start at 4 mag and are 
spaced by 4 mag. For both images, the grayscale in units of S/N. 



James et al 



2002h . The submm/mm dust emissivity is particularly im portant for star form ation 



2002t) . dust 



studies. Since molecules are known to deplete inside prestellar cores (iBergin et al. 
emission may represent the best tracer of the gas density distribution just prior to the onset 
of gravitational collapse. Thus measurements in the submm/mm define the initial conditions 
from which a core collapses to form a star Alternatively, the density distribution can be mapped 
through extinction ^ by measuring th e near-infrared colour excess towards giant stars in the back- 



ground of a cloud ( iLada et al.l 19941) . High resolution an d S/N maps can be obtained for object in 



the foreground of dense stellar fields (lAlves et all2001l) . 

This is the case for the dark cloud Barnard 68, a starless globule seen in the foreground of 



the Galactic Bulge. 



Alves et al. 



( 120011) have produced a high resolution extinction map of the 
cloud, measuring the H-K colour excess of nearly 4000 stars in its background. In this paper, we 
compare the extinction map with observations of submm/mm dust emission of similar resolution. 
Observations are described in Sect. 121 In Sect.|3lwe deri ve the dust emissivity froni the correlation 
between emission and extinction (in a way similar to iKramer et alJll998l l2002h . We will also 



adopt a temperature gradient within the cloud, which was derived from a model of dust heating. 
Finally, the derived emissivities are compared with other estimates from literature in Sect. 0] 



2. Submm/mm Observations 

Barnard 68 was observed with the Submillimetre Common User Bolometer Array (SCUBA) 



at the JCMT dHolland et al 



19991) and with the SEST Imaging Bolometer Array (SIMBA; 
Nvman et alJ2001l) . Both instruments have hexagonal arrays of 37 bolometers: SCUBA can ob- 
serve in the submm at 850 /zm with FWHM=14.5" while SIMBA operates at 1.2 mm and has a 
beamsize of 24". 
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SCUBA observed an area of 5'x5' around Barnard 68 in scan-map mode, chopping within 
the field of view. The resuhant image is convolved with the chop function, which is removed by 
means of Fourier Transform analysis . The Emerson II technique is used to minimise the noise 
I Jenness. Lightfoot. & Hollandll998 1: six scans of the field are made with diff'erent chop throws 
(20", 30" and 65" in RA and Dec). In total, we coadded eight sets of scans, six observed in 
Mai-ch/June 2002 and two in July 1998 (the latter retrieved from the SCUBA Archive). Standard 
data reduction was performed with the dedicated package SURF. After flat-fielding and masking 
of noisy bolometers, a baseline was removed from each bolometer Maps were corrected for 
atmospheric extinction (the zenith optical depth was Tg5o=0. 12,0. 15,0.3 in the 1998 run and in 
the March/June 2002 runs, respectively) and for correlated sky noise. Calibration was obtained 
from scan-maps of Uranus. The final map has a sky noise of 20 mJy beam"' (1-c), equivalent to 
3.5 MJy sr ' for the SCUBA beamsize. 

SCUBA scan-maps suffer from large scale undulations in the background of the source, 
which are due both to the p oor sampling of small spatial freq uencies and to uncertainties in 
the baseline removal (see, e.g. 



Visser. Richer. & Chandleg2002h . On our image, the undulations 



show up as a negative background around the source, of the order of the sky noise. We subtracted 
a mean background estimated on several sky apertures. The results of the next section are not 
significantly changed by the subtraction. The integrated flux at 850 //m inside the 2-cr isophote 
is F850//in=4. 0+1.0 Jy. SCUBA also took maps with the 450 jjm array. The signal at this wave- 
length is poor, because of the background fluctuations and of the larger sky opacity (T450 ~ 0.5 
for the best data). The peak emission is detected at 3-cr level (1-cr - 200 mJy beam"'; beam- 
size FWHM=10") and it is not possible to conduct the surface brightness analysis described in 
Sect.|3] Nevertheless, we could estimate the integrated flux, which is F4g()^m= 13+4 Jy. 

SI MBA observations w ere taken in fast scanning mode, without using a wobbling secondary 



20011) . The bolometer array scans the observed area with a speed of 80" 



mirror (iNvman et al. 

s~' in azimuth and a scan-to-scan separation of 8" in elevation. Spatial frequencies are converted 
into temporal frequencies and the sky signal is filtered out with a low-cut filter. Each map of 
Barnard 68 covers 600" x 312" in azimuth and elevation, and takes 7 minutes to complete. 
Barnard 68 was observed during June and October 2001. A total of 118 maps were coadded 
together, for a total integration time of 13.8 hours. The mean zenith optical depth is ti 2=0.16. 

Data were reduced with the software package MOPSI by R. Zylka. After deconvolution ac- 
cording to the frequency passband, subtracting a baseline to each scan, correcting for sky opacity 
and gain elevation the correlated sky noise fluctuations and artifacts due to electronics were re- 
moved. All maps of Barnard 68 were rebinned in a single image, with pixel size 8". In an iterative 
way, a preliminary image was used as source model to improve the removal of correlated noise. 
Uranus maps were used to calibrate the data. Night-to-night fluctuations suggest a calibration 
uncertainty of 20% (l-cr). The final SIMBA image of Barnard 68 has a residual noise of 5.5 
mJy beam"' (1-cr), equivalent to 0.36 MJy sr"' for the SIMBA resolution. The integrated flux is 

^'l.2mm=0.7±0.2Jy 
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Fig. 2. Pixel to pixel correlation between dust emission and Ay extinction, at 850 yum (left) and 
1.2 mm (right). The solid line is the isothermal fit, assuming a temperature 7'iso=12K, while 
the dashed line and the emissivity values refer to the temperature gradient model with external 
temperature rext=14K (see text for details). All high signal-to-noise points with Ay > 10 have 
been used for the fit. 



1 20011) rthe measured NIR 



Rieke&Lebofskvll985 ). 



The emission is compared to the extinction map of lAlves et al. 
colour excess has been converted to Ay using a standard reddening law; 
All images have been smoothed to the SIMBA resolution, registered together (pointing accuracy 
is better than a few arcsec for the submm/mm observations) and resampled to a common pixel 
size of 12". Fig.[Oshows the final SCUBA (left) and SIMBA (right) images, with superimposed 
Ay contours. 



3. Analysis 

The morphology of the extinction and emission maps in Fig.^is quite similar Dust emission 
with S/N>3 traces regions with extinction Ay >6 and Ay>10, in the SCUBA and SIMBA images, 
respectively. In analogy with extinction, emission comes from a round region with a south east 
tail. At 1 .2 mm, there are hints for a secondary emission peak. In Fig.|2lwe show the pixel-to-pixel 
correlation of the emission maps with Ay. As already seen in Fig.H there is a clear correlation 
between emission and extinction. For regions with Ay>10 the scatter in the correlation with the 
submm data is of the order of the sky noise (1.7 MJy sr"' for the smoothed image). For Ay<10 
the scatter increases because of the (possibly hotter) south-east tail. The pixels clearly belonging 
to the tail can be seen in Fig.|2laround Ay=8 and /v=10 MJy sr"'. The scatter at 1.2 mm is more 
uniform and slightly larger than the sky noise, because of the complex central morphology. A 
forthcoming paper will be devoted to the analysis of the central peaks and of the south-east tail. 
For optically thin radiation (as for the submm/mm emission of Barnard 68) we can write 

,,= -x^x^^LS^IU^, a, 

«, 1.086 fp dl 
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where T and p are the temperature and density of dust grains, functions of the position within the 
cloud, By is the Planck function and the integral extends along the line of sight through the dust. 
In the isothermal case, the last term on the right reduces to By{T). Otherwise, a knowledge of the 
dust temperature and density distribution is needed. 

In Fig.|3]we show the Spectral Energy Distribution (SE P) of Barnard 68. Our integrated 
fluxes at 850 //m and 1 .2 mm are similar to those presented by 



Ward-Thompson. Andre. & Kirk 



1 2002h . from which we took the ISOPHOT fluxes at llOfim and 200yum. We fitted a isothermal 
modified blackbo dy to the datapoints as suming a value for /3, the emissivity spectral index (Ky oc 
v^). For /3- 1 .5-2 ( Dunne & Bale J2OOI ), we find acceptable fits for temperatures Tiso = 11 - 1 3 K. 
We adopt Tiso = 12 + 2 K (1-cr) for the isothermal case, incorporating in the broad error the 
large uncertainties on The solid line in Fig. I^is an indicative SED for fi-l.l and risn=12 



K. Similar temperature h ave been previously obtained for Barnard 68 dWard-Thompson et al 



2QQ21 



Hotzel et al 



Using Tiso, we fitted Eq. ^to the correlations of Fig. |2] and derived 



a single cloud- averaged value for the submm and mm emissivities. We obtained /C850/im/*'v = 
3.5 + 1.0 ■ IQ-^ and kij^Jk^ = 9.0 + 3.0 ■ 10"^ Errors were estimated with a bootstrap technique 
and are dominated by the uncertainty on calibration and on temperature. 

The fitted correlations are shown as solid lines in Fig.|2l While the fit is acceptable at 1.2 
mm, a simple linear relation is unable to reproduce well the correlation observed at 850 jjm for 
Ay > 20 since the ly/Ay ratio decreases for increasing Ay. This is expected if the tempera- 
ture in the cloud core is lower than in the external part, as a result of the dust shielding of the 



external radiation field 1 Z}ig£2gii_^lri]:Sl ev. & Galli 



other clouds (see, e.g. 



Kramer et al 



20011) . Such a trend has been observed in 



1998h . A similar behaviour could also be explained with a 
isothermal dust distribution, if the dust emissivity is lower in den ser regions. However, models 
and observations su ggest that emissivities increase in dense cores ( Ossenkopf & Henningll994 ; 



Kramer et al.l2002 ) 



We derived th e temperature gradient inside the cloud from a n improved versio n of 



the Zucconi et al 



Alvesetal 



1 20011) . the 



model (Gongalves et al., in preparation). As in 
adopted density distribution is that for a Bonnor-Ebert sphere, a pressure confined isother- 
mal sphere in hydrostatic equilibriu m. Dust is heated by a local Interstellar Radiation Field 
Galli. Walmslev. & Goncalvesl2002 ). The internal dust absorption is fixed by the measured ex- 
tinction and the assumed abso rption law Ky/Ky, for which we used the tabulated values given by 



Ossenkopf & HenningI (Il994h for different models of dust coating/coagulation in dense clouds. 
Wliile the absolute value of the temperature depends on the adopted dust model (for A > 30//m, 
Ky/Ky increases going from bare grains to grains with ice coating up to grains that undergo coag- 
ulation), the temperature gradient is found to be relatively i ndependent of grai n characteristics. 



Alves et al 



20011) about 1.5 times 



with the temperature at the truncation radius Text (0.06 pc; 
higher then the core temperature. This is because the absorption law does not change signifi- 
cantly with the dust model at shorter wavelengths, where most of the absorption occurs. In the 
foUowing, we adopt the model radial gradient and derive Text from the data. 
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For the assumed density and temperature distributions, we simulated emission maps with 24" 
resolution and 12" pixels, in analogy with the data. The integrated SED is shown with a dashed 
line in Fig.|3 Again, a broad range of best fit temperatures is obtained, rext=12-15 K. We use 
Text - 14+2 K (1-cr). Finally, the simulated Ay vs ly correlation of Eqn.[2is fitted to the observed 
data to derive Ky/Ky. We obtained^ 

^850^m/^V = 4.0+1.0-10-5 
Kl.2mm/Kv = 9.0 ± 3.0 ■ 10"^ 

These correlations are also shown in Fig.|2lwith dashed lines. Although the fit to the submm data 
has improved, the emissivities are nearly the same as for the isothermal case. 



4. Discussion 



The derived sub mm/mm emissivities are shown in Fig. 0] together with estimates from 
literature. Emissivities derived from IR AS/COBE observations of diffuse Galactic dust 
Bianchi etalJl999 1 alues predicted by the popu- 



Boulanger et al 



lar 



92^ 



Ossenkopf & Henning 



Drgingj^Lgg j 12^4 1 mod el, which is shown in Fig.l^as a solid line. Dotted lines refer to the 



1 19941) models for bare grains, grains with a thin ice coating and grains 
with thin ice after coagulation has proceeded for 10^ years in a gas with density 10^ g cm"^ (from 
bot tom to top, r e specti vely) . 



James et al 



1 2002h derived the emissivity from a sample of local galaxies observed with 



SCUB A, estimating the dust mass from metals. Their value is compatible to that of lDraine & Lee 



1 19841) . indicating similar properties for diffuse dust in external galaxies and in the Milky Way. 
Our 850 um emissivity is hi gher than that of diffuse dust: a comparison with the models of 
Ossenkopf & Hennind ( 1 19941) seems to suggest that grains in Barnard 68 possess molecular ice 

^ Equivalent to Kgso^m = 1.5±0.4cm- g"' and/fi2mm = 0. 35 + 0.1 cm^ (a ssuming grain radius OA/jm, 
grain density 3 g cm"'' and V-band extinction efficiency 1.5; 



HildebranJ1983 ). 



I 




200 400 600 800 1000 

A [fj-m] 

Fig. 3. SED of Barnard 68. The solid line is the isothermal model with Tiso = 12 K. The dashed 
line is for the model with a temperature gradient, with Text = 14 K. Both cases assume /3 - 1.7. 
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Fig. 4. Submm/mm emissiviti es compared with other Uterature estimates (details in the 
text). A standard ratio 



verted to k^Ik^ as in 



K^Jjcy ( iRieke & Lebofskv 



Ossenkopf & Hennind (I1994). Data fromlDraine & Led (119841) an d 



Bianchi et al 



1 19991) and 



19851) has been used to derive k^Ikv from 



Alton et al 



James et al 



1 20021) are con- 



200(1) . Points at 850jum have been 



slightly shifted in A for ease of presentation. 



mantles or have coagulated into fluf fy aggregates. Howev er, Kssof^m/Kv is within 2cr of that for 



bare grains. In a very recent paper, 



Kramer et al 



I 2002b measured K^sOf^m/K-v for a molecular 



ridge and four embedded cores in the dark cloud ICS 146. Their radially averaged values of the 
four cores, derived assuming y6=2, range from 2.1 • 10"^ to 6.1 • 10"^ and are broadly compatible 
with our determination. At 1 .2 m m, our emissivity is si milar to that of diffuse dust but lower than 



the value previously estimated by 



Kramer et al 



( 19981) on a ICS 146 core. 



The ratio /fgso^m/'i'i.imm we derive implies a value for /3 larger than what expe cted from mea 



surements on Galactic diffuse dust and on laboratory cosmic dust analogues dMennella et al. 



19981) . However, the error on the ratio is large and /? = 2 is still compatible (within l-cr) with 
our result. Contam ination of the submm and mm fluxes by the '2CO(3-2) and '^00(2-1) lines 



I Averv et alll987^ was found to be negligible. 

Clearly, more observations are needed to narrow down the errors in the determination of 
the emissivity. The analysis presented here has to be repeated on a large sample of objects of 
relatively simple morphology like Barnard 68. Only with a statistical sample will it be possible 
to study the variation of dust properties with the environment. 

Acknowledgements. We are grateful to J. Alves, R. Cesaroni, A. Natta, C. Lada and C. &amer (the referee) 
for useful comments and discussions. 
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